Theory of Kondo lattices and its application to 
high-temperature superconductivity and pseudo-gaps in cuprate oxides 
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A theory of Kondo lattices is developed for the t-J model on a square lattice. The spin sus- 
ceptibility is described in a form consistent with a physical picture of Kondo lattices: Local spin 
fluctuations at different sites interact with each other by a bare intersite exchange interaction, 
which is mainly composed of two terms such as the superexchange interaction, which arises from 
the virtual exchange of spin-channel pair excitations of electrons across the Mott-Hubbard gap, 
and an exchange interaction arising from that of Gutzwiller's quasi-particles. The bare exchange 
interaction is enhanced by intersite spin fluctuations developed because of itself. The enhanced 
exchange interaction is responsible for the development of superconducting fluctuations as well as 
the Cooper pairing between Gutzwiller's quasi-particles. On the basis of the microscopic theory, 
we develop a phenomenological theory of Zoio-temperature superconductivity and pseudo-gaps in 
the under-doped region as well as high-temperature superconductivity in the optimal-doped region. 
Anisotropic pseudo-gaps open mainly because of d7-wave superconducting low-energy fluctuations: 
Quasi-particle spectra around (±7r/a, 0) and (0, ±vr/a), with a the lattice constant, or X points at 
the chemical potential are swept away by strong inelastic scatterings, and quasi-particles are well 
defined only around (±7r/2a, ±7r/2a) on the Fermi surface or line. As temperatures decrease in the 
vicinity of superconducting critical temperatures, pseudo-gaps become smaller and the well-defined 
region is extending toward X points. The condensation of afy-wave Cooper pairs eventually occurs 
at low enough temperatures when the pair breaking by inelastic scatterings becomes small enough. 

PACS numbers: 74.20.-z, 71.10.-w, 75.10.Lp 



I. INTRODUCTION 

It is an important issue to elucidate the mechanism of 
high-temperature (high-T c ) superconductivity occurring 
on CuC>2 planes iS* 3 ^* 5 - Cooper pairs can be bound by an 
exchange interaction^ According to an early theory^& 
one of the most possible mechanisms is the condensation 
of d'y-wave Cooper pairs bound by the superexchange 
interaction. Because the on-site Hubbard repulsion is 
strong in high-T c cuprate oxides, theoretical critical tem- 
peratures T c for conventional Bardeen-Cooper-Schricffcr 
(BCS) superconductivity or the condensation of isotropic 
s-wave Cooper pairs must be very low. As long as the su- 
perexchange interaction is antiferromagnetic and it works 
mainly between nearest neighbors, theoretical T c 's for the 
dry wave are much higher than those for other waves. 
Many experiments imply that the condensation of d'y- 
wave Cooper pairs occurs in the cuprate oxidesAiS 

The hola^ doping for which observed T c 's show a max- 
imum is called an optimal doping; dopings are classified 
into over- and under-doped ones according to whether 
they are less or more than the optimal one. Both of 
normal and superconducting states in the optimal- and 
over-doped regions can be explained within the theo- 
retical framework of the early theory^ However, super- 
conductivity in the under-doped region significantly de- 
viates from the prediction of the early theory. Ob- 
served T c 's decrease with decreasing hole dopings, al- 
though theoretical T c 's never decrease. Superconducting 
gaps £g(0) at T = K keep increasing with decreasing 
hole dopings. 12 i 13 i 14 i 15 i 16 i 17 i 18 i 19 i 20 Experimental ratios of 
6g(0)/A:bT c are much large than 4.28 given by the early 



theorytS 

Normal states above T c are anomalous in the under- 
doped region; many experiments imply or show the 
opening of pseudo-gaps in quasi-particle spectra. For 
example, the longitudinal nuclear magnetic relaxation 
(NMR) rate l/XjT becomes smaller as temperatures 
T are lowered^ According to angle resolved photoe- 
mission spectroscopy (ARPES)^2iiii^ii 5 . quasi-particle 
spectra at the chemical potential are swept away around 
(±7r/o, 0) and (0,±7r/a), with a the lattice constant. 
Scanning tunneling spectroscopy (STS) directly shows 
the opening of pseudo-gapsii 6 * 1 ^ 1 ^ 9 * 2 ^ Pseudo-gaps 
start to open around the so called mean-field critical tem- 
perature defined bjii 8 ^ 9 - T c q ~ eG(0)/4fcs, which is much 
higher than T c in the under-doped region. 

The anisotropy of pseudo-gaps deduced from ARPES 
is similar to that of e?7-wave superconducting gaps. One 
may argue that the opening of pseudo-gaps must be a 
precursor of that of superconducting gaps. On the other 
hand, the development of pseudo-gaps with decreasing 
T is not monotonic in several STS data^i 9 ^ the magni- 
tude of pseudo-gaps shows a maximum at a temperature 
between T c and T c q and it decreases with decreasing T 
below the temperature. One may also argue that pseudo- 
gaps must be suppressed by superconductivity. It is one 
of the most important issues to clarify how or whether 
the opening of pseudo-gaps above T c is related with that 
of superconducting gaps below T c . 

The cuprate oxides with no hole dopings are Mott- 
Hubbard's insulators. Even if holes are doped, they must 
lie in the strong-coupling regime defined by U/W > 1, 
with U the Hubbard repulsion and W the bandwidth. 
The Hartree-Fock approximation, the random-phase ap- 
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proximation (RPA) and their more or less improved the- 
ories such as the self-consistent renormalization (SCR) 
theory starting from RPAjSSi the fluctuation exchange 
(FLEX) approximation starting from RPA and so on are 
never any relevant approximations for the cuprate oxides; 
they are relevant only for weakly correlated electron liq- 
uids in the weak-coupling regime defined by U/W <C 1. 

Hubbard's 23 and Gutzwiller's 2 ^ are within the single- 
site approximation (SSA) and are relevant for U/W > 1. 
When we take both of them^ we can argue that the 
density of states must be of a three-peak structure, 
Gutzwiller's quasi-particle band at the chemical potential 
between the lower and upper Hubbard bands. This spec- 
ulation was confirmed in another SSA theory^ a narrow 
band, which is nothing but Gutzwiller's band, appears 
at the top of the lower Hubbard band when the electron 
density is less than half-filling. Gutzwiller's band is re- 
sponsible for metallic behavior, and the Mott-Hubbard 
splitting occurs in both metallic and insulating phases. 

The validity of the SSA theories, Hubbard's and 
Gutzwiller's, implies that local fluctuations are responsi- 
ble for the three-peak structure. Local fluctuations are 
rigorously considered in one of the best SSA's that in- 
clude all the single-site terms' Such an SSA is reduced 
to determining and solving selfconsistently the Anderson 
model>2£ which is one of the simplest effective Hamilto- 
nians for the Kondo problem. The Kondo problem has 
already been solved; 30 ! 31 ! 32 ! 33 ! 34 ! 35 ! 36 so that many useful 
results are available. One of the most essential physics in- 
volved in the Kondo problem is that a localized magnetic 
moment is quenched by local quantum spin fluctuations 
so that the ground states is a singled or a normal Fermi 
liquid. 32,33 The Kondo temperature Tk is defined as a 
temperature or energy scale of local quantum spin fluc- 
tuations. The so called Abrikosov-Suhl or Kondo peak 
in the Kondo problem corresponds to Gutzwiller's band. 
Their bandwidth is of the order of IcbTk, with fee the 
Boltzmann constant. On the basis of the mapping to 
the Kondo problem, we argue that lattice systems must 
show a metal-insulator crossover as a function of T: They 
are nondegenerate Fermi liquids atT> Tk because lo- 
cal thermal spin fluctuations are dominant, while they 
are Landau's normal Fermi liquids at T <C Tk because 
local quantum spin fluctuations are dominant and mag- 
netic moments are quenched by them. Local-moment 
magnetism occurs at T > Tk while itinerant-electron 
magnetism occurs at T <C Tk- Superconductivity can 
only occur at T <C Tk- 

The superexchange exchange interaction was originally 
derived for Mott-Hubbard's insulators! 37 ! 38 ! 39 ; 40 ! 41 One 
may suspect that it must work only in insulating phases. 
On the other hand, it was recently shown from a field the- 
oretical approach that it arises from the virtual exchange 
of spin-channel pair excitation of electrons across the 
Mott-Hubbard gap^ 2 ^ Gutzwiller's band plays no role 
in this exchange process. It is unquestionable that the su- 
perexchange interaction works even in metallic phases as 
long as the Mott-Hubbard splitting exists. An exchange 
interaction arising from the virtual exchange of spin- 



channel pair excitations of Gutzwiller's quasi-particles 
also plays a role in metallic phases^ 2 ^ 3 - 

A perturbative treatment of exchange interactions in- 
cluding the above two ones starting from an unperturbed 
state constructed in the SSA is nothing but a theory of 
Kondo lattices. Because the SSA is rigorous for Landau's 
normal Fermi liquids in infinite dimensions^ it can also 
be formulated as a 1/d expansion theory, with d the spa- 
tial dimensionality. The theory has already been applied 
to various phenomena occurring in strongly correlated 
electron liquids, not only high-T c superconductivity 4 ^ 4 ^ 
but also itinerant-electron ferromagnetism^, the Curie- 
Weiss law of itinerant-electron magnets^ field-induced 
ferromagnetism or metamagnetism, 47 itinerant-electron 
antiferromagnetism^ and so on. The early theory 7 ^ of 
high-X^ superconductivity is also within a first approxi- 
mate framework of the theory of Kondo lattices. 

Because of anomalous properties of the under-doped 
region, not a few people speculate that a novel mech- 
anism must be responsible for high-T c superconductiv- 
ity. On the other hand, we confirmed in solving the 
Kondo problem that adiabatic or analytical continuity 
is one of the most important concepts in physics^ 4 * 3 ^ 9 " 
It is, of course, one of the basic assumptions of Landau's 
Fermi-liquid theory. The theory of Kondo lattices also 
relies on this concept. There is no phase transition or no 
symmetry change between high-T c superconductivity in 
the optimal-doped region and low-T c superconductivity 
in the under-doped region. Then, analytical continuity 
tells us that the theory of Kondo lattices, which can ex- 
plain superconductivity in the over- and optimal-doped 
regions, must apply to low-T c superconductivity in the 
under-doped region. One of the purposes of this paper is 
to apply the theory of Kondo lattices to low-T c supercon- 
ductivity and pseudo-gaps in the under-doped region. 

This paper is organized as follows: In Sec. |Hj a the- 
ory of Kondo lattices is developed for the t-J model. On 
the basis of the microscopic theory developed in Sec. [HJ 
a phenomenological theory of low-T c superconductivity 
and pseudo-gaps is developed in Sec. IIHI where several 
parameters are phenomenologically determined instead 
of completing self-consistent procedures involved in the 
microscopic theory. Discussions are given in Sec. IIVI 
Conclusions are given in Sec. 

II. THEORY OF KONDO LATTICES 

We consider the t-J model on a square lattice: 
H = - ^2 *»j ffl L a jv - 5 J ' S ^ 

ija (ij) 

a^ajtaj^a^, (2-1) 

i 

with S l = i J2af3 °z & ) a i a a t0 , witn &y and 

o z the Pauli matrixes, and the summation over (ij) re- 
stricted to nearest neighbors. Here, notations are conven- 
tional. The second term is the superexchange interaction. 
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Because we are interested in the cuprate oxides, we as- 
sume that J~ — (0.10-0. 15)eV and the transfer integral 
between nearest-neighbors is t\ ~ (0.4-0.5) eV: 



J/\h 



-(0.2-0.3). 



The third term is introduced to exclude any doubly oc- 
cupied sites, so that E/oo/|fi| — > +oo. Because no phase 
transition is possible at nonzero temperatures in two di- 
mensions, we assume that weak three dimensionality ex- 
ists but the reduction of T c due to it is small. 

When Landau's normal Fermi liquid is considered, the 
single-particle self-energy is divided into single-site and 
multi-site terms: S CT (je n ,k) = E .(ig n ) + AS cr («e„,k). 
The single-site term S CT (ie„) is identical to the self- 
energy for a mapped Anderson model, which should be 
self-consistently determined and solved^ The multi-site 
term is divided into two terms: 



AX a (ie n ,k) = AS CT (k) + AS;(ie n ,k), 



(2.3) 



where the first term is independent of energies ie n and 
AE^(ie„,k) — > in the limit of |e n | — > +oo. 

First, we construct an unperturbed state in a renor- 
malized SSA, where not only T, a (ie n ) but also AE CT (k) 
are included. Then, the single-particle Green function is 
described in such a way that 
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i £„+ n - EQt) - T, a (ie n ) - AE ff (k) 



where fi is the chemical potential, and E(k) = 
—2 tii]i s (k) is the dispersion relation of electrons, with 
ti the Zth nearest-neighbor one of the transfer integral 
and rji s {k) the form factor of the Ith. nearest neighbors: 

r?is(k) = cos(k x a) + cos(kya), (2.5a) 
?72 S (k) = cos[(k x + ky)a] + cos[(fc x - k y )a], (2.5b) 

and so on, with a the lattice constant. The mapping 
condition to the Anderson model is given by2& 



^£Gg ) >(ie n ,k)=G< r (ie n ), 



(2.6) 



with N the number of lattice sites and G a (ie n ) the Green 
function of the Anderson model; the on-site interaction 
of the Anderson model should also be infinitely large. 
The single-site term is expanded in such a way that 

t„(ie n ) = Efff (l-^ 7 )te«+(l-0»)-ffAiBfl'+-.., (2.7) 

in the presence of infinitesimally small fields H, with g 
being the g factor and fin the Bohr magneton. The ex- 
pansion (|2.7|l is accurate for \s n \/kBTfc < 1; we approx- 
imately use it for \e n \/k.BTjc ^ 1- The Wilson ratio, 
which appears in the Kondo problem, is defined by 



(2.8) 



Because charge fluctuations are suppressed, W s — > 2 in 
the limit of half-filling or n — > l^O^^ with n the elec- 
tron density per site. Theories in SSA give approximate 



values of 



1/1 



24.26 



The dispersion relation of 



(2.2) quasi-particles in the unperturbed state is given by 



&(k) = ~- EQt) + t + AE CT (k) - fi 



(2.9) 



and the Green functions are given by Gi ^(z£ ra ,k) 
(l/4>y)ga (i£ n ,)&) + (incoherent part), with 



ffS 0) (fe»,k) 



1 



Co(k) +ij sgn(£„)' 



is,. 



(2.10) 



being the coherent part describing quasi-particles; the 
incoherent part describes the lower and upper Hubbard 
bands. Here, a phenomenological life-time width 7 is 
introduced; sgn(x) is defined in such a way that sgn(;r) = 
— 1 for x < while sgn(a;) = 1 for x > 0. 

When the irreducible polarization function in spin 
channels is denoted by ir s (iu)i,q), the susceptibility of 
the t-J model, which does not include the conventional 
factor \g 2 tJt B , is given by 



Xs(*^/,q) 



27r s (iw;,q) 



1 - | 5 J(q) + C/qoJ n s (iuji,q) 



, (2.4) with 



J(q) =2J ms (q). 



(2.11) 



(2.12) 



The irreducible polarization function is also divided into 
single-site and multi-site terms: ir s (iuji,q) = 7t s (iu>i) + 
An s (ioj[,q). The single-site term n s (iuJi) is identical to 
that of the Anderson model, so that the susceptibility of 
the Anderson model is given by 



2TT s (iu)l) 



1 - UooKs^LOl) 



(2.13) 



An energy scale of local quantum spin fluctuations or the 
Kondo temperature is defined by 



k B T K 



(2.14) 



Bec ause t /polxsC^) I — > +00 and f7oo|x«(M*>j,q)| 
Eq. p.llfl can also be written in another form: 



Xs(«w;,q) 



I s (iui,q)x s (iui) ' 



with 



I s (iwi,q) = J(q) +2U^ An s {iuj h q). 



(2.15) 



(2.16) 



This expression is consistent with a physical picture of 
Kondo lattices: Local spin fluctuations at different sites 
interact with each other by an intersite exchange inter- 
action. Then, we call I s (iu>i,q) an exchange interaction. 



According to the Ward relation^ the irreducible 
single-site three-point vertex function in spin channels 
is given by X s {ie n , ie n +iwi; iui) = 4> s [l-U^^Ui)] or 

\ s {ie n , ie n + iu>r,iu)i) = 2cj> s /U 00 x s {iui) 1 (2.17) 

for En/ksTK — > +0 and wi/UbTk — > +0. We approxi- 
mately use Eq. I|2.17[l for nonzero e n and w;. Then, the 
so called spin-fluctuation mediated interaction in the lon- 
gitudinal spin channel is given by 



2 _ i 

[/ooA s (0,0;0) [xs{iui,q)-x s (iuji)} = 4> S -I*(iu)i, q), 



(2.18) 



with 



1 



ll s (iut,q) 



-I*(wi,q) = - l 
4 1 - ±I s (iuji,(i)xs{toJi) 



ij 5 (iwi,q) + 



~I s (iw;,q) 



~] 2 



X«(^,q) 



(2.19b) 



In Eq. I|2.18[) . the single-site term is subtracted because 
it is considered in SSA. The interaction in the transver- 
sal channels is also given by Eq. (|2.18|1 : the spin space 
is isotropic in our system. Because of Eqs. (|2.18|l and 
(|2.19(l . we call I s (iuJi,q) a bare exchange interaction, 
J*(iw;,q) an enhanced one, and <j) s an effective three- 
point vertex function. 

The bare exchange interaction l|2.16[) is mainly com- 
posed of three terms: 

I s {iu h q) = J(q) + J Q (iu h q) -4A(^,q). (2.20) 

The first term is the superexchange interaction Ij2.12|l . 
The second term is an exchange interaction arising from 
the virtual exchange of pair excitations of quasi-particles. 
When Eq. (|2.17l) is made use of, its lowest-order term in 
intersite processes is given by 4 ^ 

1 W 2 r 

-J Q (^,q) = -^[p(^,q)-P (^)], (2.21) 



with 



P(m,q) = -fc B T^l^.g(. 0) (^n + ^,k + q) 

Ere k(T 



x 5 (°)(ie n> k). 



(2.22) 



The summation over e n can be analytically carried out 
as is shown in Eq. (| A .If) . Here, the single-site term, 
Po(iuii) = (1/-/V) J^q P{iui, q), is subtracted because it 
is considered in SSA. The third term corresponds to the 
so called mode-mode coupling term in the SCR theory. 22 
Because the w-linear imaginary term in Pq [lo + iO) is can- 
celled by that in l/x s (w + i0), 42 Eq. (|2.15|l is approxi- 
mately given by 



Xs(^ + «0,q) 



1 



1/^(0) -i/> + i0,q)' 



(2.23) 



(a.) 



(b) 



(c) 



^ I ^ H \ ? ,7<^a^2vX eJ^^jL' 



FIG. 1: Feynman diagrams for the multi-site self-energy. A 
solid lines stands for gz* (ie n ,k), a wavy line for I*(iui,q), 
and a solid circle for W s - Diagram (a) gives the sum of 
AE CT (k) and AE^"^ ( is„ , k) . All the orders of two types of di- 
agrams such as (b) and (c) are considered for AY>i sc ^ (ie n , k). 



for |o>| < ksTx', the exchange interaction is given by 
(2.19a) I' s (oj+M,q) = J(q) + J' Q (u+iO,q) - 4A(0,Q), (2.24a) 



-J' Q {uj+iO,q) 



X 2 s(0) 



[P(u+i0,q)-P o { 



(2.24b) 



with Q = (±7r/a, ±7r/a). 

As is shown in Eq. J22J), there are two types of multi- 
site self-energy corrections. One is the Fock term due to 
the superexchange interaction. When only the coherent 
part is considered, it is calculated in such a way that 



1 



^AE CT (k) = -W t 



3 ,~, 2 k B T 



N 



,fen0+„(0) 



= 2JSr] ls (k), 



ffW(ie„,k+q) 
(2.25) 



with 



3,- , knT 



4 



Vuipy^g^^p). (2.26) 



Here, the factor 3 appears because of three spin channels. 
According to Gutzwiller'sjSi <fi 7 — > +oo in the limit of 
n — ► 1 and the bandwidth of quasi-particles vanishes. 
When the Fock term is included, however, the bandwidth 
is about | J\ in the unperturbed state even in the limit of 
n — > 1, if there is no disorder. 

So far the unperturbed state is constructed in a renor- 
malized SSA where E CT (ie„) and AS (T (k) should be self- 
consistently calculated. Next, AE^.(ze n ,k) in Eq. i|2.3|) 
is perturbatively considered in terms of I s (iu>i , q) or 
I* (itoi , q) . We consider two types of corrections shown in 
Fig.[[] AS' CT (ie n ,k) = AS^ F) (ze„,k) + AZ { f c \ie n ,k). 
One arises from antiferromagnetic spin fluctuations: 



xg^(te n + iLJ h k+q), (2.27) 



with 



U AF (iwi,q) = ~\l* s {iw u q) - J(q) 



(2.28) 



Here, the superexchange interaction J(q), which gives 
the Fock term, is subtracted. 
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The other arises from G?7-wave superconducting fluc- 
tuations. We define a reducible polarization function of 
^7-wave particle-particle or Cooper-pair channel by 

U aa , (iu h q) = jr^Vdi (p) Vdj (k) / dre lu}lT 
P k ^° 

a (p+|q)^ T ) a (-p+|q)-a( T ) ffl (-k+Iq)-,' a (k+Iq),')' 

(2.29) 



with 



77rf 7 (k) = cos^a) — cos(k y a). 



(2.30) 



Cooper pairs are bound with the enhanced exchange in- 
teraction H2.19fl . It is expanded in such a way that 

q) = +2^J?(iw,)»fe(q). (2.31) 

i 

Because the nearest-neighbor plays a major role, 

we only consider it and we ignore other ones. When 
only ladder diagrams such as those shown in Figs, ^b) 
and (c) are considered, fluctuations of different waves are 
decoupled from each other, so that H aa t(iu)i, q) satisfies 
the following coupled equation: 

I1 CTCT (iuji , q) = vr d7 (iuji , q) + i vr d7 (iwj , q) 4> 2 S II CT _ CT (z , q) 

- j7r d7 (zwi,q)7*^n CTCT (iwi,q), (2.32a) 
Ilo-o-^w^q) = -7rd 7 (zwj,q)/*0fn«r(^,q) 

- i^ d7 (iw, , q)7*^n CT _ CT (to, , q) . (2.32b) 

Here, relations of n o -_ . (ia;^ , q) = n_ CT(T (za;/, q) and 
no-o-(iwi,q) = n_ cr _ cr (icji,q) as well as 

2j?i(k-p) = r?i s (k)?7i s (p) +r? d7 (k)r?d 7 (p) 

+r? x (k)r? x (p) + %(k)%(p), (2.33) 

with 7? x (k) = ^/2sm(k x a) and %(k) = v2sin(fcj,a), are 
made use of; the energy dependence of Il(iu>i) is ignored 
and J* = (Re[2J (u> + i0)}) is an average over a low-energy 
region such as |w| < k B T K ; n dy (iu;i,q) = 7r^ 7 (iw/, q)/^ 7 , 
with 

7T^,q) = -M^^ (kh (0) (fen)k+ i q) 

x^Hen-MJi.-k+lq). (2.34) 

The summation over e n can be analytically carried out 
as is shown in Eq. (|A.2Jl . As is discussed in Sec.|U T c 's for 
the d"f wave are much higher than T c 's for other waves, 
so that low-energy fluctuations of the d'j wave dominate 
over those of other waves. Then, we consider only them. 
A superconducting susceptibility for the dj wave is given 
by 



Xd-,(m,q) = 4> 2 S [n CTCT (icj;,q) -n 



r(»Wi,q)] 



l + |/*^ 2 7r* 7 (z^,q)' 



(2.35) 



The self-energy correction is given by 
^A4 SC) (^,k) = -^£^X( k +!q) 



xJ7 d7 (iw i ,q).g^ 0) (-i£ n -icj(, -k-q), (2.36) 



with 



Z7 d7 (zo;;,q) = -/j* xLC^.q). 



(2.37) 



In Eqs. H2.35|) and l|2.37[) . the factor 3 appears together 
with I* because of three spin channels. 

The Green function renormalized by antiferromag- 
netic and superconducting fluctuations is given by 
G CT («£ n ,k) = g a (isni k)/0 7 + (incoherent part), with 



g a (ie n ,\s) 



' - J-AS;(ie n ,k) 



-i -l 



^ 0) (ze„,k) ^ 



(2.38) 



with AE' a {ie n ,k) = AE& 4F) (ie„,k) + AEi 5C) (ie„,k). 
The density of states for quasi-particles is given by 



with 



Pk( £ ) = Im 



g a (e + iO, k). 



(2.39) 



(2.40) 



III. SEMI-PHENOMENOLOGICAL THEORY 

A. High-r c superconductivity 

On the basis of the formulation in Sec. [D] we develop 
a phenomenological theory. First, we consider the un- 
perturbed state. The dispersion relation £o(k) of quasi- 
particles in the unperturbed state is expanded as 



£o(k) = M* 



E 



2i*rfo(k). 



(3.1) 



Although the expansion coefficients such 

and so on should be self-consistently determined, we treat 

them as phenomenological parameters. We assume that 



t\ > 0, t% = -0.3^, 



(3.2) 



and we ignore other tVs. Eq. (|3.2|l is consistent with ex- 
periment. Furthermore, we treat t\ as an energy unit for 
the sake of simplicity although it depends on tempera- 
tures, electron densities and so on. 

The density of states for quasi-particles in the unper- 
turbed state is given by 



(3.3) 
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According to the Fermi-surface sum rule^i the electron 
density is given by that of quasi-particles, so that in the 
limit of T -> +0 and 7^+0 

» = 2^E eiE " 0+ ^ 0) (^,k) 

- 2jdep 1=0 {e)f 1 {e) = 2fde Pl {e)f- f=0 {e), (3.4) 



with 



/ 7 ( £ ) = ^ + - Im 



1 2 27rfc B T 



(3.5) 



with the di-gamma function. Note that p 1= o(e) = 
(l/N) E k S[e - Co(k)] and / 7=0 ( e ) = l/[e^ T + 1]. 

We assume Eq. even for nonzero T and 7. 

Nonzero 7's are introduced, partly, for the convenience 
of numerical processes. Fig. Ufa) shows p 1 ie) for several 
7. As long as 7's are rather small such as 7/|i*| < 0.3, 
we expect that essential features are never wiped out by 
such 7. Fig. H£b) shows n as a function of /it*, Fig.[2Ic) 
shows p* as a function of n, and Fig. Hfd) shows Fermi 
surfaces for 7 = and several n. 

According to Fermi-liquid relationship and the map- 
ping condition l|2.6|) . it follows that [xs(0)] T ^ +0 = 

2W s pj(0) in the limit of 7 — > 0. In almost half-filling 
cases (n ~ 1), W s ~ 2, as is discussed in Sec. [H] Ac- 
cording to Fig. EJa), it follows that p 7 (Q) ~ 0.2/|i*| 
and ksTx — 1.2|i*|, unless the chemical potential is in 
the vicinity of the two-dimensional van Hove singularity 
or 7's are extremely large. The specific heat coefficient 
due to local spin and charge fluctuations is given by^P 
7c = (2/3)/c|<V 7 (0). 

We take the following phenomcnological expressions 
for susceptibilities: 



Xs(u + i0, Q + q) 



Xd 7 (w + i0,q) - , , 



n 2 + (qa) 2 - iuj/T s 



n d + (qa) + c u uj - iu/T, 



(3.6a) 
(3.6b) 



with Q = (±7r/a, ±7r/a). 

According to an analysis of experimental data^ 
Xs(0, Q) — 75 states/eV, 1/k s ~ 2.3 and uj$f — F s k 2 ~ 
14 meV; n s = 0.43, T s = 74 mcV and X s(0, Q)T s n 2 = 
1.1. In the optimal-doped region, the specific- heat coef- 
ficient is about 14 mJ/mol K 2 , so that \t*\ ~ 50 meV. 
Then, we assume 



Xs (o,Q)r sK 2 



1, 



r s /|tti~i, (3.7) 



for the optimal-doped region. 

When the dispersion relation of quasi-particles is given, 
it is straightforward to calculate polarization functions. 
For example, Fig. shows P(uu + iO, q). Then, we can 
estimate some of parameters appearing in Eq. (|3.6a() . Ac- 
cording to Eqs. (|2.23|) and (|2.24a|l . it follows that 



Xs (0,Q)F sK 2 
1 



w 2 
x 2 M 



Im 



dP(uj + i0,0) 



-J 



W 2 
X 2 (0) 



duj 

9P(0,Q + q) 



d{qaf 



j=0 



q=0 



(3.8a) 
(3.8b) 



Parameters estimated from Fig.[3]together with Eq. I|3.8(l 
for 7/|^| = 0.3, n = 0.9 and k B T/\t\\ = 0.1 are consis- 
tent with those given by Eq. (|3.7|l . 

Fig. 0] shows 7T^ (to + iO, q). There exists an w-linear 
real term in l/ir^ioj + iQ, q), as is shown in Fig. Efc). 
We ignore it and we assume that c w = in Eq. (|3.6b() : 
nonzero c w can never play any crucial role in physi- 
cal properties examined in this paper. According to 
Eqs. (jOS) and l|3~6b")l . it follows that 



W 2 



xS 7 (o,o)r d7K 2 7 
w 2 



= -Im 



d 



d 



1 



9w7r| 7 (o;+i0,0) 



, (3.9a) 



j=0 



(3.9b) 



q=0 



x5 7 (o,oK 7 L 9 M 27r S 7 (o,q)_ 

Then, it follows from Figs. 0^c) and (d) that 

X * dl (0,0)T dl K 2 dl ^W 2 /2, T dl /\tl\ ~ 1/2, (3.10a) 
for 7/|iH = 0.3, 7i = 0.9 and fc B T/|i*| = 0.1; 

x2 7 (o,o)r d7K 2 7 c^ 2 /3, r d7 /|^| ~ 1/20, (3.10b) 

for 7/|^| = 3,n= 0.9 and k B T/\t\\ = 0.1. 

Note that there are crucial differences between effects 
of 7 on spin and Cooper-pair channel fluctuations: J 
appears in Eq. I|3.8|l but it does not in Eq. I|3.9|l : the u>- 
linear imaginary term becomes small with increasing 7 in 
the spin channel, as is shown in Fig. 01a), while it does 
not in the Cooper-pair channel, as is shown Fig.^fc). For 
larger 7's, r s 's are larger but F smaller. Small 

r^'s play a crucial role in the formation of pseudo-gaps 
studied in Sec. ImBl 

When we take Eq. i|3.6|) . the energy integration in 
Eqs l|2.27l) and (|2.36|) can be analytically carried out: 
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FIG. 2: Single-particle properties of the unperturbed state for T — K, t* > and = — 0.3t*. (a) Density of states for 
quasi-particles p T (s), (b) the electron density n as a function of the effective chemical potential /i*, and (c) /x* as functions of 
n. In these three figures, solid, dotted, broken, dot-broken and two-dot-broken lines are for y/\t*\ = 1CT 3 , 0.1, 0.3, 1 and 3, 
respectively, (d) Fermi surfaces for 7 = and 19 electron densities such as n = 0.1 x i, with 1 < i < 19 being an integer. Dotted 
lines show those for n —0.5, 1.0 and 1.5. More than half-filling cases (n > 1) can be treated, if we take the hola— picture. 




-] -Of ii Chi I ftj. L M 



FIG. 3: Spin-channel polarization functions: (a) P{u + «0, Q)|t*| and (b) P(+i0, Q + q)|t* [, with Q = (±7r/2, ±7r/2), for 
n — 0.9 and fcsT/|t*| = 0.1; 7/|t*| = 0.3, 1 and 3 are assumed, respectively, in (1), (2) and (3). In Fig. Eta), solid and dotted 
lines show the real and imaginary parts, respectively; In Fig.^b), solid and dotted lines show those for q || (1, 0) and q || (1, 1), 
respectively; the polarization functions are almost isotropic around Q. (c) [P(+i0, q) — Po(+iO)] \t* | for k B T/\tl | = j/\t* | = 0.1, 
and (d) [P(+i0,q) — Po(+*0)] |i* | for k B T/\tl\ = 7/|£i| = 0.3; solid, dotted, broken, dot-broken and two-dot-broken lines are 
for n — 1, 0.9, 0.8, 0.7 and 0.6, respectively. Because P(+i0,q)'s have peaks at wave-numbers different from Q for small n, 
effective k s 's had better been estimated from Figs.[3Jc) and (d) instead of Fig. |3Jb). 



with 



J-AS[ jtJ?) (e + i0,k) 



^AI]^ C) (e + i0,k) 





1 


= 9AF^ 




1 


= gsc^ 




"1 


- 9sc 


4 



q|<9c 



|q|<?c 



1 5( E ,-e(k + q),r d7 [^ 7 + ( ga ) 2 ] 



|q|<?c 



9af = 3 



1 2 



-I s (uj + i0,Q)W s 



(3.11a) 
(3.11b) 
(3.11c) 



(3.12) 



S(e,x,V) 



k B T 



+ 



1 

+ 2^ 
1 

+ 2^ 



T e + 27 - x + iT 2tv 



2nk B T 2 



-ie + T 1 
2-Kk B T + 2 



2irk B T 



2irk B T 



—is + ix + 7 
2nk B T 



£ + 17 — x — iT 



2nk B T ' 2 



— ix + 7 1 
2nk B T + 2 



1\ , , /— ie+ix+y 
"' + ^( 2,k B T 
1 



£+17 — x + iT 



£ — 17 — x + ir 



(3.13) 
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FIG. 4: Polarization functions of d-y-wave Cooper-pair channel: (a) 7r2 7 (w + i0, 0)|i*|, (b) 7r^ 7 (-|-i0, q) |tj |, (c) l/7Td 7 (o' + iO, 0)|f* | 
and (d) l/7r2 7 (-M0, q) |tj for n = 0.9 and ksT /\t\\ = 0.1. In Figs.|lja) and (c), solid and dotted lines show real and imaginary 
parts, respectively. In Figs. |4jb) and (d), solid and dotted lines show those for q || (1,0) and q || (1,1), respectively; the 
polarization functions are almost isotropic. In (1), (2) and (3), 7/|£*| = 0.3, 1 and 3 are assumed, respectively. As 7's increase, 
the w-linear imaginary term of l/7r2 7 (oj + iO, 0) increases while its g-quadratic term decreases. 

Because only low-energy antiferromagnetic fluctuations are crucial in Eq. (|3.11a|l , we use instead of Eq. I|2.28|) 



3 r 1 3 

l?AF(a;+i0,q) = - iJicj^q) - J(q) + — I^(iuji, q)x s («w;, q) ~ 3 

4 L J 4^ 



4^.(0, Q) 



X s (w+i0,q). 



(3.14) 



Because we use the phenomenological forms, we restrict 
the summation over q in Eqs. I|3.11|l and we assume that 



q c a — 7r/3. 



(3.15) 



Because q c is rather small, we approximately use 
Eq. (|3.11c|l instead of Eq. lj3.Hb(l . In this approximation, 
there are two interesting properties: When k s = and 
T s = r</ 7 , it follows that 

— AE^)(+i0,k x ) = -— \A^ c \ +l 0,k x )Y , 
9af gsc L J 

(3.16a) 

for X points defined by kx = (±7r/a, ±n/a). In deriving 
Eq. (|3.16a|) . two relations of 5(0, x, T) = S*(0,-x,T) and 
£(k x + q x ,k y + q y ) = £(fc x + q y ± 7r/a, k y + q x ± n/a) for 
k = kx are made use of. The other is 

AE^°)(£ + i0,k) = 0, (3.16b) 

for k || kfcf, with kj\/ defined by k A / = (±7r/a,0) and 
(0, ±7r/a). It is certain that the renormalization by dr/- 
wave superconducting fluctuations is anisotropic. 

The electron density is crucial for the nesting of the 
Fermi surface and the spin susceptibility. Because we use 
the phenomenological form H3.6afl for the spin suscepti- 
bility, however, the electron density itself is never crucial. 
Then, we assume n — 0.9 here and in the following part of 
this paper. Fig.EJa) shows AE^ F) (e+iO, kx), Fig. Eft)) 
shows AEi Ai " ) (e+iO,k A f/2), and Fig.^c) and Fig.EJshow 
AT,^f G \e+i0, kx). Fig.dshows the imaginary parts of 



their static components (e = 
widths of quasi-particles. 
Because |t*| ~ 50 meV, 



-i0) or inelastic life-time 



in the optimal-doped region. If only the superexchangc 
interaction as large as Eq. (|2.2[1 is considered, it follows 
that |J a (0,Q|)/4|tJ| > \J\/\t\\ = (2-3) and \Ij\ /\tt\ > 
\J\/\t\\ = (2-3). When we take Eqs. and ftj.lOaL 

Qaf = (12-27) VF 2 |^| 2 and g sc = (4.5-10)1^ | 2 . If 
we take W s ~ 2, g AF ~ (48-100) |^| 2 and g SC - 
(72-160)|t|| 2 . On the other hand, observed temperature 
dependence of resistivity implies 

(l/0 7 )lm AZ a (+i0, k) ~ -(l-2)* fl T, (3.18) 

in the optimal-doped region, so that gAF and gsc should 
be much smaller than these. For example, if we take 



W s (0.7-1), 



(3.19) 



gAF — 8(t|) 2 and gsc — 2 (if) 2 . Then, it follows from 
Figs. Ha) and (b) that 

(l/0 7 )lm AEi AF) (+i0, k) ~ -k B T, (3.20a) 

for k = kx and k = kjf/2; it follows from Fig.[7{a) that 

(l/^ 7 )Im AY,^ (+»0, kx) =s -k B T. (3.20b) 

Eq. H3.20|l is consistent with Eq. I|3.18|l . As long as W s ~ 
2, we cannot explain observed resistivity of the optimal- 
doping region. 

There is a similar drawback for theoretical super- 
conducting critical temperatures T c . According to 
Eq. I|2.35|l . T c 's are given by 



l + -7fW s V(0,0) = 0. 



(3.21) 
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(3.17) 



The imaginary part of the self-energy has a reduction 
effect of T c . Because its static part plays the most crucial 
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FIG. 5: Self-energy corrections: (a) {is + iO, kx)\t*\/<f>ygAF and (b) T,i AF \ie + iO,]i.M /2)\ti\/4>jgAF for n = 0.9, 

K S = tt/8 and T s /\tt\ = 1; (c) E<- SC) (is + iO, kx)|t* \/4>y9sc for n = 0.9, k Ai = tt/8 and r d7 /|^j = 1. In (1), (2), (3) and (4), 
ksT j\t\ \ = 0.1, 0.3, 1 and 3 are assumed, respectively; solid and dotted lines show the real and imaginary parts. 




FIG. 6: Self-energy corrections: Tr„ C \ie + iO, kx)|i* \/4>ygsc for n = 0.9, Kd-y = tt/8 and smaller r d7 's; (a) F d7 /|ti| = 1/4, 
(b) r dj /\tt\ = 1/8, and (c) r d7 /|t*| = 1/16. In (1), (2), (3) and (4), k B T/\tl\ = 0.1, 0.3, 1 and 3 are assumed, respectively; 
solid and dotted lines show the real and imaginary parts. 



role, we consider it as the pair breaking. Because it is 
almost linear in T as is shown in Fig. \7\ we assume a 
phenomenological form: 



7 = c o 



Cl 



(l- Cl )-r4(k) 



k B T. 



(3.22) 



Fig. |H1 shows T c for three cases of c\ = 0, 1/2 and 1 
as a function of cq. The reduction of T c is as large as 
T c o/T c ~ 2 for Co — (1-2), with T c q critical temperatures 
in the absence of any pair breaking. If we take 



I* W s 



1*11 = 



-(1.5-3), 



(3.23) 



observed T c in the optimal-doped region can be ex- 
plained. Eq. (|3.23|l is satisfied when we take Eq. (|3.19|) . 

In the formulation in Sec. ^ g„ (i£ n ,k) appear 
in Feynman diagrams. These (is n , k) had better 
be renormalized into g a (ie n ,]s.) in an improved theory. 



When this rcnormalization is made 



45 



is renormalized 



into </> 7 (k) and W s into W s (k) = </> s /0 7 (k). The energy 
derivative of the intersite self-energy is estimated from 
results shown in Figs. |S] and |SJ so that Wg(k) must be 
significantly smaller than W s ~ 2. 53 In this paper, how- 
ever, we simply treat W a as another phenomenological 
parameter; we assume Eq. (|3.19|) . that is, W s ~ (0.7-1). 



B. Pseudo-gaps and low-T c superconductivity 



The following four properties are crucial for the 
under-doped region, (i) Observed superconducting gaps 
£g(0) at T = OK increase with decreasing hole 
dopingsi 12 ' 13 i 14 ' 1 ^ 1 ^' 1 Ti 1 ?i 1 ?i 2l i l It implies an increase of the 
Cooper-pair interaction, that is, the enhanced exchange 
interaction If. It is also theoretically argued in Sec. II VI 
that II increases with decreasing hole dopings. (ii) As 
is discussed in Sec. IIII Al the energy scale of supercon- 
ducting fluctuations T e i-y is smaller for a larger life-time 
width 7 of quasi-particles. (iii) The life-time width 7 is 
larger for smaller Td- ( or when low-energy superconduct- 
ing fluctuations are developed, as is shown in Fig. (iv) 
As hole dopings decrease, cuprate oxides are closer to an 
antifcrromagnetic phase and inelastic scatterings by anti- 
ferromagnetic fluctuations more substantially contribute 
to 7. Then, we claim that the under-doped region can be 
characterized by large 7 and small Td-y, or large life-time 
widths of quasi-particles and low-energy superconducting 
fluctuations. 

An increase of 7 causes a decrease of r c j 7 and the de- 
crease of Td-y causes an increase of 7. Because of this co- 
operative effect, an increase of /* and a development of 
antifcrromagnetic fluctuations with decreasing hole dop- 
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FIG. 7: Inelastic life-time widths: Figs. 0a) and (c) show Im[AEi' 4F ' (+i0, kx)] |i* \/4>y9AF, which is equal to 
— Im[AE' sc '' (+i0, kx)l *\ti\/4>~igsc when their corresponding parameters are the same between the two channels; Figs. 0b) 
and (d) show Im[AE^ AF '(-|-iO, kjw"/2)] \t\\/^gAF- In each figure, n = 0.9 and 7/|t*| = 0.3 are assumed; solid, dotted, broken, 
dot-broken and two-dot-broken lines are for r s /[i*| = 1, 1/2, 1/4, 1/8 and 1/16, respectively. In Figs. 0a) and (b), k s is 
assumed to be independent of T: (1) k s — 7r/4, (2) n s — 7r/8, and (3) k s = n/12. In Figs. (c) and (d), re s is given by 
Eq. (1) K = tt/4, (2) k = tt/8, and (3) k = tt/12. 
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FIG. 8: Superconducting T c as a function of the T-linear co- 
efficient c : (1) 7 = cofc B T, (2) 7 = c fc s T [l + r?i,(k)/4] /2, 
and (3) 7 = cofcsT?7j 7 (k)/4. Solid, dotted, and broken lines 
show T c for I* Ws /\t{ \ = — 1, —2, and —3, respectively. 



ings causes an increase of 7 and it eventually causes a 
large decrease of Td-y and a large increase of 7. Al- 
though this cooperative effect should be examined in a 
self-consistent way, we treat Td^ as another phcnomeno- 
logical parameter; we assume that it is as small as 



1/64 < W|^| < 1/2. 



(3.24) 



in the under-doped region. 

The life-time width 7 makes the us and q dependences 
of P(us + i0,q) small, as is shown in Figs. Eta) and (b). 
Then, the us dependence of l/x s (w + «0,q) is small in 
the under-doped region. However, its q dependence can 
never be as small as that of 1 /x2 7 (w + *0, q), because the 
q dependence of the superexchange interaction is never 
suppressed by 7. Therefore, r s 's must be large rather 
than small in the under-doped region. However, we as- 
sume for the sake of simplicity that r s /|t*| = 1 even in 
the under-doped region. 

Although gAF and gsc f° r the under-doped region are 
larger than those for the optimal-doped one, we assume 



similar values to those for the optimal-doped one: 

g AF = 8(tl) 2 , gsc = 2(^) 2 . (3.25) 

Figs. Eta)-(d) show relative intensities of quasi- 
particles at the chemical potential, p£(0), in a fourth of 
the Brillouin zone for four cases such as (a) r^/jt* | = 
1/64, (b) lV|t;| = 1/32, (c) = 1/16, and 

(d) r d7 /|^| = 1/2; other parameter are 7 = 0.3|t*|, 
k B T/\t\\ = 0.1 and k s = K dy = n/8. When T dl are 
small, spectral intensities are swept away and are almost 
vanishing around X points, (±7r/a, 0) and (0, ±7r/o), as is 
shown in Fig. Eta). The vanishing of quasi-particle spec- 
tra around X points is because of large life-time widths. 
When r^ 7 are large, on the other hand, life-time widths 
are small and the Fermi surface or the Fermi line forms 
a closed line, as is shown in Fig. Eld). 

Fig. HOf a) shows the density of states for quasi- 
particles, p*(e), together with p 7= o.3|tj| (e) for compar- 
ison, as a function of e for four values of r<2 7 : r,j 7 /|f*| = 
1/64, 1/32, 1/16, and 1/2; other parameter are the same 
as those for Fig. El Pseudo-gap develop as r<2 7 decrease. 

Figs. 1111 and 1121 show relative intensities of p^(e) as a 
function of e and k. Parameters used for Figs.llllandll2l 
are the same as those used for Figs. Eld) and (a), respec- 
tively; Fig. 1111 presumably corresponds to the optimal- 
doped region and Fig. 1121 to the under-doped region. A 
weakly dispersive band around X points in Fig. Illf a) 
is because of saddle points in the dispersion relation of 
quasi-particles. The spectra shown in Figs. lllf aWd) are 
rather normal because life-time widths are small. On the 
other hand, we can see interesting and anomalous prop- 
erties in Figs. Et a )-(d) because life-time widths large. 
Large pseudo-gaps open around X points, and almost 
dispersionless peaks or flat bands appear along T-X be- 
low and above the chemical potential, as is shown in 
Fig-Eta). Such flat bands much below the chemical po- 
tential must correspond to observed flat bands^ Pseudo- 
gaps become smaller as wavenumbers are closer to kjif/2, 
as is shown in Figs. I12f b)-(d). 
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FIG. 9: Relative intensities pk(0) as functions of k x and k v for n = 0.9, kBT/\t\\ = 0.1, k s — Kd-y = (w/8), r s /|£* | = 1, 
g AF = %\tl? and g SC = 2\tt\ 2 : (a) r^/ltf] = 1/64, (b) T d7 /|^| = 1/32, (c) r*y/|tl| = 1/16, and (d) T d7 /j^| = 1/2. 
When Td'y are small, quasi-particles are well denned only around k = (ir /2a, it /2a). Lighter parts have large intensities. The 
well-defined region extends toward to (w/a, 0) and (0, 7r/o) as Fd 7 becomes larger. 



(a) (b) 




FIG. 10: Density of states p*(e). In Fig. HUf al. parame- 
ters are the same as those for Fig. [5] solid, dotted, broken 
and dot-broken lines correspond to Fig. Q^a), (b), (c) and 
(d), respectively. In Fig. HOf b). we assume that r s /|t*| = 1, 
r d7 /|*!| = 1/64 and n s = K dj = {ir/8)y/Wk^T/\q\; solid, 
dotted, broken and dot-broken lines are for kBT/\t*\ = 0.3, 
0.2, 0.1 and 0.05, respectively. In both of Figs. HoT a'l and (b), 
p 7 (e) for 7/|t*| = 0.3 is shown in a two-dot-broken line for 
comparison. In Fig. HOf a). pseudo-gaps become larger with 
decreasing Fd 7 . In Fig. HOf bl. they become smaller with de- 
creasing T. 



IV. DISCUSSIONS 

In an improved theory, the effective transfer integral 
t\ of quasi-particles is further renormalized by antiferro- 
magnetic and superconducting fluctuations. In this pa- 
per, however, we treat t\ for the unperturbed state as one 
of phenomenological parameters and we take it as the 
energy unit. Therefore, T-dependences discussed in this 
paper are qualitative. 

Low-T c superconductivity in the under-doped region 
can be explained by large inelastic life-time widths of 
quasi-particles. For example, (l/</> 7 )Im AS CT (+i0, kx) — 
— (10-20)fcsT for parameters used to obtain Fig. As 
is shown in Fig. |SJ the reduction of T c is very large for 
such large life-time widths. 

When superconducting gaps as large as ec{T) open, 
spin and superconducting low-energy fluctuations with 
M ^ € g{T) are depressed and the pair breaking caused 
by the fluctuations must also be depressed. There- 



fore, the development of £g(T) must be very rapid 
with decreasing T. This argument is consistent with 
experiment It also implies that the reduction of 
ec(0) by inelastic scatterings must be very small, so that 

£G(0) = ££(0)^0 

kBT c ksTco T c T c ' 

with T c o critical temperatures in the absence of any pair 
breaking. In Sec. IIII Al we argue that T c o/T c ~ 2 in the 
optimal region, so that 6g(0)/A:bT c ~ 8. This number 
8 is consistent with experiment^ It is quite reasonable 
that e G (0)/k B T c > 8 or e G (0)/k B T c > 8 in the under- 
doped region. 

It is likely that n 2 s and n 2 dl decrease almost linearly in 
T with decreasing T and they vanish at critical tempera- 
tures, Tjy and T c , for antiferromagnetism and supercon- 
ductivity respectively. Then, we assume that both of T/v 
and T c are zero or very low, so that 

k s = K ^10k B T/\tl\, K dl = n^lQk B T/\tl\. (4.2) 

We ignore the T-dependence of T s and T^. Figs. Etc) 
and (d) show inelastic life-time widths of quasi-particles 
for three cases of kq = 7r/4, 7t/8 and 7r/12. It decreases 
sublincarly in T rather than linearly. It is likely that 
resistivity varies in T a with a < 1 in the under-doped 
region. It is desirable to develop a microscopic theory in 
order to predict a precise number of the exponent. 

One may argue that the opening of pseudo-gaps is evi- 
dence for the formation of preformed Cooper pairs above 
T c . If this scenario is the case, pseudo-gaps must increase 
with decreasing T at any temperature region. Exper- 
imentally, however, pseudo-gaps decrease with decreas- 
ing T at a temperature range close to T c m^^ This ob- 
servation contradicts the above scenario. Instead, it is 
a piece of evidence that the opening of pseudo-gaps is 
because spectral intensities of quasi-particles around X 
points are swept away by inelastic scatterings. Fig- ITUT b) 
shows p*(s), together with p 7 =o.3|t5;|( e ) f° r comparison, 
as a function of e for four cases such as k B T/\t\\ = 0.3, 
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FIG. 11: Relative intensities of Pk(e)- Parameters are the same as those for Fig. |SJd), and they correspond to those for the 
optimal-doped region except for the electron density n = 0.9. Lighter parts have large intensities. 
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FIG. 12: Relative intensities of Pk(e). Parameters are the same as those for Fig. Eta), and they correspond to those for the 
under-doped doping region; the electron density is assumed to be n = 0.9. Lighter parts have large intensities. 



0.2, 0.1, and 0.05; we also assume Eq. (I4.2f) with k = n/8 
as well as T 8 /\tl\ = 1 and T d7 /|i*| = 1/64. The decrease 
of Kd-y causes an increase of pseudo-gaps, if other param- 
eters are constant. In Fig- HclT b). however, the magnitude 
of pseudo-gaps decreases with decreasing T because the 
decrease of life-time widths dominates that of Kd~/- This 
result is consistent with experiment 

It is reasonable that K^L is much smaller below T c q 
£g(0)/(4/cb) or that the coefficient 
is quite different between 
below and above T c0 . If this is actually the case, pseudo- 
gaps start to open around T c q. It is desirable to study 
the T-dependences of Tdj and Kd-y in order to predict the 
T-dependence of pseudo-gaps. 

When Eq. l|3.6a|) is assumed, the longitudinal NMR 
rate is given by 



than it is above T c q 
of the linear T term in K dl 



— — oc — lim Im 



Xs(^ + »0,q) = Xs(0,Q) 

lu r.. 



Because life-time widths are large in the under-doped 
region, r s 's are also large there. When the opening of 



pseudo-gaps is mainly caused by superconducting fluc- 
tuations, Xs(0, Q) must be reduced and r s must be en- 
hanced. It is interesting to examine if the development 
of superconducting fluctuations can actually explain ob- 
served reduction of 1/TiTiSl It is difficult to explain the 
reduction of Xs(0,q) an d \/T\T by the development of 
antiferromagnetic spin fluctuations. 

Experimentally, the so called kink structure is observed 
in the dispersion relation of quasi-particlcs determined 
by ARPES experiment^ Very tiny kink structures can 
be seen around the chemical potential in Figs. Hir e) 
and I12f c): they are too tiny to explain ARPES exper- 
iment. When pseudo-gaps or superconducting gaps are 
developed, spectra of antiferromagnetic and supercon- 
ducting fluctuations may have certain gap- like structures, 
which are often called resonance modes. It is interest- 
ing to examine if such resonance modes actually exist in 
(4.3) Xs{u + i0, q) and Xd 7 ( w + *0, q). One may argue that 
what cause pseudo-gap must cause the kink structure. 
It is interesting to examine effects of superconducting 
fluctuations extending in the wavenumber space, whose 
energies are rather high, as large as \zbTk or of the order 
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FIG. 13: (a) Nearest neighbor Re[|t*|Pi(w + iO)] and (b) 
next-nearest neighbor Re[|t*|p2(oJ + iO)] as a function of u 

-- 1. 




FIG. 14: (a) Nearest neighbor Re[Xi(w)] and (b) next- 
nearest neighbor Re[X2(w)]. Note that the directions of k s 
and UJ axes are opposite to each other between (a) and (b). 



of the bandwidth of quasi-particles. 

Cooper pairs of the dj wave are mainly bound by the 
enhanced exchange interaction. It is divided into the bare 
one and the enhanced part, as is shown in Eq. I(2.19b(l . 
As is shown in Eq. 1(2.20(1 . the bare one is further divided 
into three terms. The first term is the superexchange 
interaction. It is antifcrromagnetic, and it is as large as 



|j/t;i ~ (2-3). 



(4.4) 



The second term is the exchange interaction arising from 
the virtual exchange of pair excitations of quasi-particles: 



1 2W 2 + °° 

-J q (oj + iO, q) = ]T P,(^ + <0)»h.(q), (4.5) 



with 



Pi{u + iO) 



N ^ 
q 



P(u + iG,q)rfo(q). 



(4.6) 



Fig. I13f a) shows nearest- neighbor Re [Pi (oj + zO)] , and 
Fig- H3f b) shows next-nearest-neighbor Re[P2(a-> + z'O)] . 
The nearest-neighbor one is antiferromagnetic for almost 
half-filling and is stronger as the electron density is closer 
to half-filling. However, it is less effective than the su- 
perexchange interaction: 



2W 2 |Pi(+iO)| 
X 2 S (0) l*J| 



< 0.1. 



(4.7) 



It is interesting that Re [Pi (uj + iQ)l is ferromagnetic for 
n < 0.7. The ferromagnetic exchange interaction must 
play a role in the reduction of T c in the over-doped re- 
gion. The third term is the so called mode-mode coupling 
term — A(u> + iO, q). This term works as an repulsive in- 
teraction for c?7-wave Cooper pairs. The enhanced part 
of Eq. 1(2.1 9 b(l is approximately given by 



[|^(0,Q)] 2 X s (^ + zO,q), 



(4.8) 



with 



,(w + iO,q) = x.(0,Q)/^5^Jf < (« + iO)» fe (cO, (4.9) 



Vis(<C 



N ^ k 2 + (qa) 2 - iuj/T s 



(4.10) 



Fig. 1141 shows nearest-neighbor Re \X\ (u> + iO)] and next- 
nearest- neighbor Re[X 2 (o; + iO)] . The nearest- neighbor 
one is also antiferromagnetic. As the electron density 
is closer to half- filling, k s 's become smaller so that the 
enhanced part is stronger. It is also less effective than 
the superexchange interaction: 



4^(0, Q) 



• x s (Q,QK 



|Re[Xi(w+iO)]| < 0.5, (4.11) 



for \lu\ ~ k B T K . We can argue from Eqs. 1)4.4(1. 14. 7|) 
and ((4.11(1 that the superexchange interaction is the main 
part of the pairing interaction in high-T c cuprate oxides, 
and that the total pairing interaction is larger as the 
electron density is closer to half-filling. 

In the weak-coupling Hubbard model, the spin- 
fluctuation mediated pairing interaction is given by 



U 2 Xs (LO + iO,q). 



(4.12) 



Not a few people claim that Cooper pairs must be 
bound by Eq. 14.12(1 in high-T c cuprate oxides. However, 
this pairing mechanism cannot apply to the cuprate ox- 
ides, which certainly lie in the strong-coupling regime. 
Eq. ((4.12(1 relevant in the weak-coupling regime is cer- 
tainly smoothly connected with the pairing mechanism 
by Eq. 1(2.18(1 or 1(2.19(1 relevant in the strong-coupling 
regime. However, Eq. 1(4.12(1 is physically different from 
Eqs. ((2.18|) and ((4.8(1 . even if it is similar to them in ap- 
pearance. Because Eq. ((4.1211 arises from the virtual ex- 
change of low-energy pair excitations of quasi-particles, 
it can also be called an exchange interaction. If we call 
the superexchange interaction a spin-fluctuation medi- 
ated interaction, on the other hand, it sounds impertinent 
because the energy scale of spin fluctuations responsible 
for the superexchange interactions is as large as the Hub- 
bard repulsion U and is much larger than the effective 
Fermi energy of quasi-particles. 



V. CONCLUSION 

A theory of Kondo lattices is developed for the t- 
J model. On the basis of the microscopic theory, a 
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phenomenological theory of high-T c superconductivity 
in cuprate oxides is developed; parameters are phe- 
nomenologically determined instead of completing the 
self-consistent procedures involved in the microscopic 
theory. 

In the strong-coupling regime, electrons are mainly 
renormalized by local quantum spin fluctuations so 
that the density of states is of a three-peak structure, 
Gutzwiller's quasi-particle band at the chemical potential 
between the lower and upper Hubbard bands. Two bare 
exchange interactions between quasi-particles are rele- 
vant: One is the superexchange interaction, which arises 
from the virtual exchange of spin-channel pair excita- 
tions of electrons between the lower and upper Hubbard 
bands, and the other is the exchange interaction aris- 
ing from that of spin-channel pair excitations of quasi- 
particles themselves. Gutzwiller's quasi-particles are fur- 
ther renormalized by antiferromagnetic and supercon- 
ducting fluctuations caused by intersite exchange inter- 
actions. The sum of bare exchange interactions including 
the two are enhanced by low-energy antiferromagnetic 
fluctuations into the enhanced one. The condensation of 
c?7-wave Cooper pairs bound by the enhanced exchange 
interaction is responsible for high-T c superconductivity 
in the optimal-doped region. 

In the under-doped region, inelastic scatterings by 



d'y-w&ve superconducting low-energy fluctuations cause 
large life-time widths of quasi-particles around X points; 
they are almost linear or sublinear in T. Anisotropic 
pseudo-gaps open because spectral intensities of quasi- 
particles around X points are swept away by strong in- 
elastic scatterings. The development of pseudo-gaps with 
decreasing T cannot be monotonic; their magnitude must 
decrease at temperatures close to T c . Superconductiv- 
ity eventually occurs when the pair breaking by inelas- 
tic scatterings becomes small enough at low enough T"s. 
Superconducting gaps ec(0) at T — K are never much 
reduced, because low-energy antiferromagnetic and su- 
perconducting fluctuations are substantially suppressed 
by the opening of ec(0) 's themselves and, in particular, 
because their thermal fluctuations vanish at T = K. 
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APPENDIX 



The summations over e n in Eqs. (|2.22l) and (|2.34l) can be analytically carried out. When the analytical continuation 
such as ilui — ► lu + iO is made, they are given by 



p(u> + *o, q) = -L R ^ £( k + q)> £( k )]> 



(A.l) 



1 



tt2> + iO, q) = — J2 *4(k)i* [w, £ (k + |q) , -£ (-k + iq)] , 



(A.2) 



with 



w - £i + 6 



1 2 2irk B T 



1 2 2nk B T 



-n- 2 + 



w - fi + 6 



-iuj - i& + 7 
2Trk B T 



i2j 
—icj — «£ 2 



1 2 2nk B T 
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^ + *£i + 7 
2itk B T 



-i(,2 + 7 
2ixk B T 



2nk B T 



\2 2irk B T ) 



(A.3) 
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